Computer Simulation of Cytoskeleton-Induced Blebbing in Lipid Membranes * 
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Blebs are balloon-shaped membrane protrusions that form during many physiological processes. 
Using computer simulation of a particle-based model for self-assembled lipid bilayers coupled to an 
elastic meshwork, we investigated the phase behavior and kinetics of blebbing. We found that blebs 
form for large values of the ratio between the areas of the bilayer and the cytoskeleton. We also 
found that blebbing can be induced when the cytoskeleton is subject to a localized ablation or a 
uniform compression. The results obtained are qualitatively in agreement with the experimental 
evidence and the model opens up the possibility to study the kinetics of bleb formation in detail. 



I. INTRODUCTION 

It is well established that morphological deformations 
of the plasma membrane may be induced by the interplay 
between the cortical cytoskeleton (CSK) and its lipid bi- 
layer (LB) [1112]. Many cells exhibit transient exoplasmic 
protrusions, known as blebs, during several physiologi- 
cal processes, including cytokinesis, cell motility, apopto- 
sis, and virus uptake pF-'F . Recent experiments revealed 
that blebs can be induced in suspended fibroblasts when 
these are subjected to a localized ablation of the corti- 
cal CSK fTUl . Blebbing is also observed when tension on 
the cortical CSK is increased after activation of myosin- 
II motors |10[ fTT] . This letter presents a computational 
study of membrane blebbing using a particle-based model 
for a self-assembled LB coupled to an explicit CSK. 

Blebs are believed to be triggered by a biochemical in- 
hibition of myosin-II leading to either a localized rup- 
ture of the cortical CSK or its detachment from bi- 
layer |10l [T^l - fTF| . They then grow into spherical pro- 
trusions, up to 2 /im in diameter, which are devoid of 
acto-myosin meshwork [15 . In non-apoptotic cells, blebs 
retract due to recruitment of actin-membrane linkers to 
the bleb's membrane, polymerization and contraction of 
actin within the bleb |T3]. Red boold cells (RBCs) were 
also shown to exhibit blebbing during their aging. In- 
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deed, a large fraction of an RBCs plasma membrane 
is shed into small vesicles that are devoid of spectrin 
CSK ^HH]. The small size of RBC-shed vesicles sug- 
gests that the precursor blebs may have a size about the 
corral size of the RBC spectrin meshwork [15 . Blebbing 
in RBCs is believed to result from a reduction of ATP 
levels during the aging process leading to a contraction 
of the sepectrin meshwork |ll8j. 

Despite its importance to various biological processes, 
few theoretical studies were performed to understand 
blebbing [TOl [T9l [20] . A theoretical or computational 
study of cell blebbing is obviously very difficult due to 
the vast complexity of cells. Simplified models are thus 
needed to infer some aspects of blebbing. Here, we pro- 
pose a particle-based model system, where the only in- 
gredients taken into account are the LB and a simplified 
CSK. Our results clearly demonstrate that membrane 
blebbing results from a subtle interplay between contrac- 
tile tension of the CSK meshwork and the LB elasticity. 

II. MODEL AND METHOD 

Our study uses a mesoscale implicit-solvent model, 
recently developed by us [2T], of self-assembled lipid 
molecules with soft interactions. Here, lipid molecules are 
coarse-grained into semi-flexible amphiphiles composed 
of one hydrophilic particle (h) connected to a chain of 
two hydrophobic (t) particles. The self-assembled one- 
component lipid vesicles do not possess a spontaneous 
curvature. The inner side of the vesicle is underlined 
with a semi-flexible polymer meshwork tessellated by tri- 
angles formed by linking vertices with semi-flexible poly- 
mer chains. The CSK configuration is composed of 162 
vertices, corresponding to A/'cor = 320 triangular corrals. 
150 vertices have six links while 12 have five links. Each 
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link is composed of 8 to 22 monomers. All links are iden- 
tical for a given vesicle. The vertices are anchored to the 
membrane through bola lipids mimicking the anchoring 
protein complexes in RBCs. We note that the CSK topol- 
ogy remains conserved during a simulation. A pictorial 
presentation of the model is shown in Fig. [l] 




Figure 1: (Color online) Lipids particles are represented by 
dots. The cytoskeleton particles are yellow (lightest) beads. 
The anchoring (bola) lipids have green (middle gray) head 
particles and blue (darkest gray) hydrophobic beads. 



The interaction potential between particles includes 
two-body interactions, harmonic interactions for the 
bonds within a lipid or CSK meshwork, and three-body 
interactions to account for bending rigidity of lipid par- 
ticles and CSK [H] 



+ E^bcnd(i-.-i,r„r,+i), (1) 



where Fj is the position of particle i, rij = jr^ — r^l, and 
Ui = h, t, or c for a lipid head, lipid tail, or a CSK 
bead, respectively. In Eq. ([l]), [/^'"^ is ^ soft two-body 
interaction between two particles i and j and is given by 
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The self-assembly of the lipid chains in this model is 
achieved through the addition of an attractive interac- 
tion between tail particles. We choose U^^^ = if either 
a or /3 = /i, and U^f^ < if a /? = t. U^^/^^ > for any 
values of a or /3. The interaction between the cytoskele- 
ton particles and lipid particles is assumed repulsive. 

The potential J7bond ensures the connectivity between 
two consecutive monomers in a lipid chain or cytsoskele- 
ton and is given by 

Kon,{r)='^ir~a,)\ (3) 

where k^^nd is the bond stiffness coefficient and is the 
preferred bond length. U^^^^^ is a three-body interaction 
potential ensuring bending stiffness of the lipid chains 
and is given by 

TTa (r- y r- \ — ^bcnd ( na ' \ 

2 V ri^^^ir.i^i+i J 

(4) 



where k^^^^^ is the bending stiffness coefficient of a lipid 
chain or CSK and 6'^ is the preferred splay angle and is 
chosen as 180° for all triplets. The values of the interac- 
tion parameters are given by 
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and rc — 2r„i, and — O.lrm- 

The thickness of the lipid bilayer in the present model 
is given 4.1r„j. The thickness of a in the fluid phase and 
in vitro condition is about 4.5 nm. We thus estimate 
r,„ ~ 0.91 nm. The area per lipid for a tensionless bi- 
layer in the present model is about 0.65 r^. The area 
per lipid for a dipalmitoylphosphatidylcholine (DPPC) 
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bilayer in the fluid phase is about 0.63nm^ fF?. There- 
fore, we estimate ~ Inm. Phospholipid bilayers in 
the fluid phase and in vitro typically have a diffusion co- 
eflicient D « 1 - 10 x lO'^^ m'^/s [13]. In our model, 
D « 0.015 r'^/r. Therefore, the time scale r « 1 — 10 ns. 
Many of our simulations were performed over durations 
of milliseconds. We note that the thickness of the lipid 
bilayer of the plasma membrane in vivo is not known, 
although we do not expect it to be much different from 
that in vitro. We also note that in vivo, the diffusivity is 
about one order of magnitude smaller. 

The particles are moved using molecular dynamics 
with a Langevin thermostat [3T]: j^i{t) = Vi(t) and 
mVi{t) = —ViU — T-Vi{t) + Wi(t), where m is the mass 
of a particle (same for all particles). F is a bead's fric- 
tion coefficient, and Wi(t) is a random force originat- 
ing from the heat bath, and satisfies (Wi(<)) = 0, and 
(W, {t) • {t')) = 6fcBrr<5, j 5 {t - t'). Equations of mo- 
tion are integrated using the velocity Verlet algorithm 
with r = ^/6m/T where the timescale r = rm{m/ey^'^ 
with rrn and e being the length and energy scales. Simula- 
tions are performed at k^T — 3e with At = 0.02r. Many 
simulations were run up to 8 x 10^ time steps, which cor- 
responds to time scales around 10ms. A large number of 
vesicles were investigated with number of lipids ranging 
between 35 000 and 2.5 x 10^, corresponding to vesicles' 
diameter between 70 and 160 nm. 



III. RESULTS 

We first focus on the phase diagram of the system (cf. 
Fig. [2|, which is described in terms of the rest area of a 

CSK corral, W^^^^ (defined as the corral area with min- 
imum elastic energy), and the area mismatch parame- 
ter, s = ^lip/^^yjQ, where ^lip is the area of the LB 
normalized by the number of corrals, A/'cor- Fig- [2] in- 
dicates a phase transition line, denoted by s*{A^{^), 
between a phase where the vesicle is conformationally 
homogeneous and a phase where the vesicle is blebbed. 
Within the homogeneous phase, the CSK conforms to 
the vesicle, while in the blebbed phase, the bleb is de- 
void of the CSK. To understand this phase behavior and 
the monotonic decrease of s* with .A^y|.Q we use an ear- 
lier argument by Sens and Gov [T^]. In the homoge- 
neous phase, the free energy of the vesicle (assumed to 
be spherical) is dominated by the curvature energy of 
the bilayer and the elastic energy of the CSK patches, 

J-h ~ StTK + (l/2)AAcore4yto (Ayto/4yto " l) . whcrC 

K is the LB bending modulus, e is the corral stretch 
modulus and y^cyto is the of a single stretched 

corral. However, if the vesicle exhibits a single bleb, 
the CSK is stress-free, and the vesicle's free energy 
is then dominated by its curvature energy, J-]^ ~ 

IGttk (l - A'fyl^/SirRl - A^^l^/SnRfj , where i?„ and Rt 

are the radii of the main part of the vesicle (with the sub- 



jacent CSK) and the bleb, respectively, both assumed to 
have almost spherical shape. In the expression of energy 
of the blebbed vesicle, it is assumed that the area of the 
bleb's neck is that of a relaxed corral. The transition 
line, in the case of a vesicle with a large number of cor- 
rals (A/'cor ^ 1), obtained by equating the free energies of 
the homogeneous and the blebbed phases, is then given 

by s* « 1 -I- IGTTK/A'core^cyto- Therefore, the transi- 
tion line decays monotonically with y^cyto- This analysis 

also predicts that s* ~ l/l^^^, where ll^to is the linear 
size of a stress-free coral. This seems to be in agreement 
with our results shown in the inset of Fig. [2] although 
simulations on much larger vesicles is needed to verify 
whether this linear relation will hold still. 
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Figure 2: (Color online) Blebbing phase diagram for closed 
vesicles in terms of corral rest area, .4^yto, and mismatch pa- 
rameter, s. The inset shows the transition line, s* , versus 
1/^cyto, where l^yto is the linear size of a stress- free CSK cor- 
ral. The solid line is a fit to the equation for S* given in the 
text. A^yjjj is varied by changing the number of beads per 
link in the cytoskeleton. 

To further investigate the nature of blebbing, the CSK 
elastic energy is shown vs. the area mismatch, s, in Fig. [3] 
for the case of A, 



(0) _ 

cyto 



92r,„, together with snapshots of 



the CSK for different values of s. This figure indicates 
that the CSK potential energy is minimized at s^un ~ 1 
slightly larger than 1. Smi,-, is slightly above 1 due to the 
slight curvature of the LB. When s < Smin, the CSK is 
compressed due to the small area of the vesicle, while 
for Smin < s < s* , the CSK is stretched, implying that it 
exerts a compressive stress on the LB as s is increased to- 
wards the phase transition. In this regime, the increased 
elastic energy of the CSK is compensated by the low cur- 
vature energy of the vesicle which is almost spherical. 
For s > s*, the CSK elastic energy becomes too high 
if it were to conform to the large membrane. Instead, 
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Figure 3: (Color online) CSK energy, per strand, vs. mis- 
match parameter, s, for the case of .4^y].o — 92r^ Inset shows 
the coral area, normalized by its rest area, vs. s, for same 
system. Blue (most right) arrow next to (d) points to the 
location of the bleb's neck. 



the free energy is minimized when the vesicle protrudes 
a bleb that is devoid of CSK. The CSK then conforms 
to an effective vesicle with smaller area. In this case, the 
CSK adopts a stress-free configuration (snapshot (d) of 
Fig. [3| which is very similar to the case of a non-blebbed 
vesicle with a stress-free CSK (snapshot (b) of Fig. 3|. 
The sharp discontinuity of the CSK elastic energy implies 
that the transition is first order, in accord with Ref. |19j . 

We now turn to the kinetics of bleb formation and 
growth following a localized ablation of the CSK. The nu- 
merical experiment shown here was performed on an ini- 
tially homogeneous vesicle with corral rest area A^^^^ — 
43. Ir^ and mismatch parameter s = 1.90 (right below 
the transition line). At t = 0, an arbitrary vertex is dis- 
sociated from its six links, and the kinetics of the system 
is then monitored over a long period of time. A snapshot 
time series is shown in Fig. [4] (see as well Movie 1 in Sup- 
plementary Material), together with a graph depicting 
the speed of the bleb's apex, Ubicb = dAh/dt {Ah is the 
distance between the bleb's apex and its base), the diam- 
eter of the bleb's neck and the average length of a CSK 
strand. The snapshot series indicates that the membrane 
caps immediately after, and at the same location where, 
the ablation is made. The cap grows rapidly during early 
times, as shown by the Wbiob- The speed reaches a max- 
imum, then decays during later times. The qualitative 
behavior of the apex's speed is similar to that reported 
by Charras et al. [12^. The diameter of the neck evolves 
non-monotonically with time: The sudden localized dis- 
sociation of one vertex from its six neighboring strands 
causes an imbalance of forces on the six neighboring ver- 
tices leading to their displacement outward and then an 
increase of the size of the strands in the ablation region. 



Figure 4: (Color online) Top panel: Snapshot of a vesicle un- 
dergoing blebbing as a result of a localized ablation of the 
CSK. The bleb occurs exactly at the location where the dam- 
age is made. Bottom left panel: Top to bottom snapshots 
of the vesicle taken from inside, at f = , 1500, and 7500r, 
respectively. The arrow at f = indicates the location of the 
ablated vertex. The dark red (gray) hallow corresponds to 
the bleb's neck. Only lipid head particles and CSK particles 
are shown. Bottom right panel: Bleb front velocity in red 
(bottom curve), the average length of a CSK strand between 
two vertices in blue (middle curve), and the average length of 
a CSK strand in the bleb's neck region in green (top curve). 



Meanwhile, this is accompanied by a compression of the 
overall CSK, demonstrated by a decrease in the average 
length of CSK strands. After the initial brief expansion 
of the neck, during which the bleb has a shape of a cap, 
the neck starts to retract. The onset of retraction of the 
CSK strands in the neck region coincides with the tran- 
sition of the bleb's shape to that of a bud. Interestingly, 
this transition roughly coincides with a slowing down in 
the growth rate of the bleb's apex. We associate this 
slowing down to the fact that the net current of lipids 
flowing into the bleb is proportional to the perimeter of 
the neck. Subsequent kinetics is characterized by slower 
growth of the bleb and a slow relaxation of the CSK's 
elastic stress. Blebs growth ceases once the CSK reaches 
its stress-free state. It is noted that no blebbing is ob- 
served following ablation if the initial value of s is much 
below the transition line. We note that the blebs in our 
study are much smaller than those reported experimen- 
tally [12 j. A quantitative comparison between the speeds 
from our study and that from experiment is premature 
since small blebs cannot be resolved in experiments. 

We also made a series of computer experiments where 
a second ablation is subsequently applied at a location 
diametrically opposite from the first ablation. In Fig. [s] 
the height of the second bleb, A/12, is plotted vs. time, 
together with the height of the growing bleb, Ahi, when 
a single ablation is applied. When two ablations are 
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Figure 5: (Color online) Left panel: Height of the second bleb 
when a vesicle subjected to two ablations. Top curve (black) 
in graph corresponds to the case of a a single ablation. Curves 
from top to bottom, excluding the first one, correspond to 
a second ablation at At = 250t, 1500t, 2000r, 2500r, and 
5000r, respectively. Right panel: A snapshot time series of 
the vesicle for the case of At = 2000t. Snapshots from top to 
bottom are taken at 7500, 15000, 50000 and 160000t, respec- 
tively. 



applied simultaneously, then two almost identical blebs 
form simultaneously and grow at the same rate (data 
not shown). However, when the second ablation is ap- 
plied subsequently, different kinetics of the second bleb 
is observed depending on the time lapse, Ai, between the 
two ablations. If At is very short (^ 250r), the second 
bleb forms and grows (cf. red curve in Fig. [5|. However, 
when At > 1500r, a second bleb forms and grows dur- 
ing intermediate times (green curve and snapshot series 
of Fig. |5|. During later times, the second bleb then re- 
tracts and eventually completely disappears (see Movie 2 
in Supplementary Material for the case of At = 2000r). 
The life time of the second bleb is reduced as the time 
interval At is increased. We note that energetically, one 
bleb is more favorable than two blebs. Therefore, even 
for small At, the second bleb should be metastable. No 
discernible second bleb is observed when the second ab- 
lation is applied after At ^ 5000r (cf. purple line in 
Fig. [5]). These results further substantiate that the ob- 
served blebbing is the result of the interplay between the 
CSK elasticity and the lipid bilayer elasticity. We em- 
phasize that in cells, the hydrostatic pressure plays an 
important role on blebbing, which is not accounted for 
in the present model. Therefore, the explicit variation 
of blebbing in our study, in terms of the parameter s is 
expected to be different from that in cells jTU] . 

We also looked at a vesicle undergoing a sudden uni- 
form contraction of the CSK, starting from an equilib- 
rium homogeneous state. This simulates the activity of 
myosin motors in generating a contraction of the actin 



Figure 6: (Color online) Top left panel: A vesicle, composed 
of 4 X 10'' lipids undergoing blebbing during a uniform contrac- 
tion of CSK to a point in the phase diagram corresponding 
to .4cyto ~ 112rJ^ and s = 2.25. Top right panel: Number of 
blebs vs. time. Inset shows same data in a double logarithmic 
plot. Bottom panel: Sequence illustrating coalescence of two 
blebs (indicated by arrows on top left panel). 

meshwork. A configuration of the vesicle during blebbing 
is shown in Fig. [6j Right after the sudden CSK's con- 
traction, many blebs are formed throughout the vesicle. 
Fig. [6] indicates a rapid increase in the number of formed 
blebs during early times, which then reaches a maximum, 
and eventually decays at later times. The late time kinet- 
ics proceeds through coalescence of blebs and retraction 
of some blebs, albeit most of blebs coalesce. An exam- 
ple illustrating the coalescence of two neighboring blebs 
is shown in the snapshot series in Fig. [6] (also see Movie 
3 in Supplementary Material). Interestingly, the number 
of blebs is found to decay as A'biobs ~ with a ~ —0.44 
which is close to 1/2. Note that a large number of runs 
and longer simulations are needed to unambiguously ex- 
tract this exponent. Since the system is undergoing spin- 
odal decomposition between the blebbed and homoge- 
neous regions. The kinetics of blebs' coalescence is sim- 
ilar to that of a two-component membranes undergoing 
phase separation without hydrodynamics for which the 
exponent of the number of buds, a = —1/2 [24]. 

IV. CONCLUSIONS 

In conclusion, we presented a numerical study of bleb- 
bing based on a model of a self-assembled LB with ex- 
plicit CSK. In this model the solvent is accounted for 
implicitly, and therefore volume constraint and cytosol 
flow are ignored here. These effects have recently been 
put forward as essential ingredients of blebbing. From 
investigation of the phase behavior of blebbing, we found 
that for small mismatch parameter, s , the equilibrium 
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state is that of a homogeneous vesicle with the CSK con- 
forming to the entire vesicle. However, for relatively large 
s, the equilibrium state is that of a blebbed vesicle, with 
the bleb being devoid of the CSK. The transition value 
of s decreases linearly with the linear size of the CSK 
corral, in agreement with the recent theory of Sens and 
Gov [19' . Blebs are observed when the membrane is sub- 
jected to a localized disruption of the CSK or a uniform 
contraction, in line with experimental observations. 
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